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ABSTRACT 


This paper encompasses the mathematical derivations of the analytic and generalized formula and recurrence 
relations to find out the radii of n umber of circles inscribed or packed in the plane region bounded by circular 
arcs (including sectors, semi and quarter circles) & the straight lines. The values of radii obtained using analytic 
formula and recurrence relations have been verified by comparing with those obtained using MATLAB codes. 
The methods used in this paper for packing circles are deterministic unlike heuristic strategies and optimization 
techniques. The analytic formulae derived for plane packing of tangent circles can be generalized and used for 
packing of spheres in 3D space and packing of circles on the spherical surface which is analogous to distribution 
of non-point charges. The packing density of identical circles, externally tangent to each other, the most densely 
packed on the regular hexagonal and the infinite planes have been formulated and analysed. This study paves the 
way for mathematically solving the problems of dense packing of circles in 2D containers, the packing of spheres 
in the voids (tetrahedral and octahedral) and finding the planar density on crystallographic plane. 
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1. Introduction 


In general, the circle packing is the arrangement of circles of equal or unequal radii on a surface (plane or curved, 
finite or infinite) without any overlap. The circle packing problems here are concerned with the packing of n 
number of external tangent circles in the plane region bounded by the circular arcs and the straight lines inside the 
square of known side, the sector of given radius and central angle and the regular hexagon of minimum size. The 
heuristic algorithms were used to pack the identical circles inside the plane container in form of triangle, square, 
rectangle, semicircle, and circle to maximize the radius for given size of container and to minimize the size of 
container for given radius [2]. The generic global optimization software packages were used to solve various circle 
packing problems and their industrial applications were reviewed [3]. The methods used here for circle packing 
are all deterministic and formulated analytically to compute the exact size of nth packed circle. The maximum 
density of the identical circles packed on the plane surface is obtained by arranging the centres of the circles in 
the hexagonal fashion which has been mathematically proved by theorem of Thue [4]. The most densely packed 
identical circles of a given radius can be enclosed by a regular hexagon of minimum side which can be determined 
in terms of number of circles and their radius. 


Similarly, the three external tangent circles of given radii a,b, and c can be enclosed by a unique circle of 
minimum radius R which has been analytically determined [1]. 


Conversely, if the radius of circumscribing circle R and the radii say a, and b are known then the radius c of third 
external tangent circle can be determined by solving the above equation for radius c in terms of a,b, and R. 
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1.1. Derivation of the radius of inscribed circle (one of three external tangent circles) 


We know that the radius (R) of circle which internally touches three external tangent circles with radii say a,b & c 
(as shown in figure-1) is given by the following generalized formula [1] 


= abc 
2/abc(a+b+c)-—(ab+bc+ca) 


Where a,b &c (€ => b =>c > 0) always satisfy the following inequality 


(0 < a,b,c < R) 


ss ab 
(va + vb) 


€ 


Now, suppose that the radius R of circumscribed circle with centre O and 
the radii of any two inscribed circles say a & b with centres A & B 
respectively are given or known then the radius c of third inscribed circle 
with centre C (see figure-1) can be easily computed by solving above 
generalized equation for unknown radius c as follows 


7 abc 
- 2: /abc(a +b+c)— (ab + be + ca) Figure-1: Three externally touching circles 


with radii a, b & c are inscribed by the circle 
with center O and R 


2RVabc(a+b+c) =abc + R(ab + bc + ca) 
Taking the squares of both the positive sides of above equation, 
2 2 
(2R abc(a+b+ c)) = (abc + R(ab+be+ ca)) 
4R?2ab(a + b)c + 4R*abc? = (aR + bR + ab)*c? + a*b?R? + 2abR(aR + bR + ab)c 


((aR + bR + ab)? — 4abR?)c? + (2abR(aR + bR + ab) — 4abR?(a + b))c + a*b?R* = 0 


((aR + bR + ab)? — 4abR?)c? — 2abR(aR + bR — ab)c + a*b*R? = 0 


Solving above quadratic equation for c using quadratic formula as follows 


2abR(aR + bR — ab) + [(2abr (aR + bR — ab)) — 4((aR + bR + ab)? — 4abR2)a2b?R2 


a 2((aR + bR + ab)? — 4abR?) 
2abR (ar +bR —ab +.,/(aR + bR — ab)? — ((aR + bR + ab)? — 4abR*)) 
~ 2((aR + bR + ab)? — 4abR2) 
_ abR(aR + bR — ab + V4abR? — 4a?bR — 4ab?R) 
= (aR + bR + ab)? — 4abR2 
abR (ar + bR—ab+2,/abR(R —a =) 
ea a = 


(aR + bR + ab)? — 4abR? 


In above expression it is worth noticing that aR + bR — ab > 0 & (aR + bR + ab)? — 4abR? > 0. The radius 
c is a positive real value (i.e. can’t be an imaginary number) therefore (R — a — b) = 0 . Thus there are two 
possible positive real values of unknown radius c which depends on the following two cases 


Case-1:R-—a—b>0OorR>a+b [discriminant is positive] 


There are two positive distinct real values of unknown radius c i.e. there 
are two different inscribed circles of radii c, & cz for given values of 
a, b & R (as shown in figure-2) given as follows 


abR (aR + bR - ab +2 abR(R — a— b)) X 
(aR + bR + ab)? — 4abR? > 


Figure-2: Two different inscribed circles 
with radii c,& cz are possible if R > a+b 

The circle with radius c, is 
inscribed/packed in the smaller region 
bounded by circles with radii a,b and R . 


Cz, C2 = 


Case-2:R-—a-—b=O0OorR=a+b (discriminant is zero] 


There are two positive equal real values of unknown radius c i.e. there are two identical inscribed circles of 
equal radii c, & c, for given values of a,b & R. In this case, the centres A 
and B will lie on the diameter of outer/circumscribed circle with centre O 
(as shown in figure-3). The equal values of radius c are given as follows 


abR(aR + bR — ab) 


“v2 = (aR + bR + ab)? — 4abR? 


> Cy, =C2 


Figure-3: Two identical inscribed 
circles with radii c,& cz are possible if 
R=a+b 


Case-3:R-—a-—b<0OorR<a+b [discriminant is negative] 
In this case, there are no real values of c, & Cz therefore this case is discarded. 


Now, it’s worth noticing that the value of radius c of inscribed circle with centre C will be less than given radii 
a,b & R (see figure-3 or 4 in above two possible cases) only when we take the negative sign in the equation 
given by above case-1 (-- -ve sign always gives decreasing value of radius c of packed circle) i.e. 


abR (aR + bR- ab-2 abR(R — a—b)) 


Ca = (GR + BR 4 ab)? — 4abR? VY Cmin < min(a, b, R) eisoass CL) 


The above formula i.e. Eq(1) should be used for computing the minimum value of radius c of (smallest) circle 
with the centre C which is inscribed/packed in the smaller region bounded by the tangent circles of given radii 
a,b & R (as shown in the above figure-2). The above eq(1) gives the decreasing value of c which is applicable in 
all such cases of circle packing as discussed below. 


2. Packing of circles in the plane bounded by square, semi and quarter circles 


Consider a square ABCD of each side x which has a semi-circle with diameter AB and a quarter circle with centre 
at the vertex B (see figure-4). 


D 


Now, consider the circles with centres C,,C,(C3,...... C, and radii 
ie a eee T, tespectively snugly fitted or packed in the plane region 
bounded by the square ABCD, semi-circle and the quarter circle (as shown 
in the figure-4). Let’s first find the radius 1, of first inscribed (blue) circle 
then the radii of all other inscribed circles 12,15, ...... T, can be easily 
computed by using mathematical relation i.e. Eq(1) as derived above. We 
have to find out radius of nth inscribed circle (labelled with blue colour). 


Now, join the centre C, of first blue inscribed circle to the centre E of semi- Figure-4: Infinite no. of small 
circle and to the vertex B of square ABCD. Drop the perpendiculars C,G and blue circles can be packed or 

i ‘ : inscribed in the plane region 
C,F from the centre C; to the sides BC and AB respectively (as shown in the pounded by square, semi and 
figure-5 below) quarter circles. 


In right AC,GB, We have, C,G=7, and C,B = x —71, (see figure-5). D 
Now using Pythagorean Theorem as follows 


(CB)? = (C,G)* + (GB)? 


=> GB=./(C,B)2 — (C,G)? 


GB = V(x —77)?-(%)? = Jx? —2xr, 


In right AC; FE, We have (see figure-5) 


x Figure-5: The blue circle of radius r, 
= ABS xe — axr , C,E= 2 +o and center C, is the first largest circle 
inscribed in the region bounded by 
square ABCD, semi-circle with center 
E and quarter circle with center B. 


FE = EB — FB = EB— C(,G 


Now using Pythagorean Theorem in right AC, FE as follows 


(C,E)? = (CF)? + (FE)? 


Gn) = (=a) + Gn) 
x2 2 


2 ae - 2 
a +xr1 =x ical er giclee — xT, 


x 
x?-4xr,=0 > abe Vx #0 


m= sbaitatten (2) 


Above Eq(2) gives the radius of first largest inscribed circle (labelled with blue colour in the above figure-5) in 
terms of edge length x of square ABCD. 


Coordinates of the centre (C,): The coordinates of the centre C, of first inscribed circle w.r.t. the vertex A 
(assumed to be the origin) of the square ABCD (see above figure-5) can be determined as follows 


x 3x 
x — coordinate = AF = AB— FB=x-7r, =x -—~=— 


4 4 


X Xx 
y — coordinate = C,F = GB = /x? — 2x7, = |x? — 2x: a. va 


; es : 3 
Therefore, the coordinates of the centre C, of first inscribed circle w.r.t. the vertex A are (=, =) (see the figure- 


5 above). Similarly, the coordinates of other inscribed circles can be computed after computing their radii. 


Radius of second inscribed circle (72): Now, draw the complete circles by extending the semi-circle with centre 
E and quarter-circle with centre B (as shown in the figure-6 below). It is worth noticing that three externally 
touching circles with centres C, ,C, and E touch the circumscribing circle (i.e. extended form of quarter circle in 


the above figure-5) with centre B. This case is similar to the case analysed in the above section 1.1 (see the above 
figure-1). 


Figure-6: The semi-circle and the quarter circle are extended to draw 
complete circles with centers E & B and radii x/2 & x respectively. 


Now, the radius r, of the blue inscribed circle with centre C, touching the circles with centres C,, E & B and the 
radii 71, = x/4, x/2 & x respectively, is found out by substituting a= 7, =x/4, b=x/2 and R = x in the 
above Eq(1) (i.e. minimum value Cyin) as follows 


abR (ar + bR — ab —2,abR(R—a =) 
(aR + bR + ab)? — 4abR? 


Cmin = 


i (5 | x5 D 
8 \-8 8 # (6292). 
= : = & __ — (5 — 2v2) 
Ge an 17x 17 
(Fr) oo 64 


Above Eq(3) gives the radius of second inscribed circle (labelled with blue colour in the above figure-6) in terms 
of edge length x of square ABCD. 


Radius of third inscribed circle (13): Similarly, consider three externally touching circles with centres 
C, , C3 and E which internally touch the largest circle with centre B (as shown in the figure-6 above). This case is 
again similar to the case analysed in above section 1.1. Thus, the value of radius rz of blue inscribed circle with 
centre C3 depends on the radii 72,x/2 & x of circles with centres C,, E & B respectively. Thus, the radius 73 of 
circle with centre C3 can be computed by substituting a = 7, b = x/2 and R = x in the above Eq(1) as follows 


abR (ar + bR — ab — 2,fabR(R —a— b)) 
min = (aR + bR + ab)? — 4abR? 


>13 
iy(8- 202)-§-x (Fy (S- 200) 45-2 35 (5 -2V) 5-2] (5-2V0)-§-x (2-5 - 21-9) 
(35 (5 — 2v2)-x+2-x +B (5 —2V2)-3) -4- (6-22) 3 x2 
x (33 - 20V2 + 85 — 34V2 — 2(5 — 2V2)V'85 — 34V2 — 66 + 402) 


289 + 297 — 180V2 — 170 + 68V2 


x (118 - 54V2 — 2(5 — 2v2)(3v2 + 1)) 


16(26 — 7V2) 
_ x(33-20V2)  x(33-20V¥2)(26+7V2) x(578-289V2)  x(2- 2) 
~ 4(26—7V2)  4(26—7V2)(26+7V2) 4678) 8 
r3 = =(2 — V2) = 0.0732233047033631xX tate (4) 


Above Eq(4) gives the radius of third inscribed circle (labelled with blue colour in the above figure-6) in terms of 
edge length x of square ABCD. 


Radius of forth inscribed circle (74): Similarly, the radius 7, can be computed by substituting a = 173, b = 
x/2 and R = x in the above Eq(1) as follows 


abR (ar + bR —ab —2,fabR(R —a— b)) 
ay = (aR + bR + ab)? — 4abR? 


>" 


§O-VD-F (FOV) x45 x-R(2-v)-$-2 5@-vD)-¥-x(x-F(2- v)-3)) 
(E(2-v2)-x+%-x+8(2-v2)-3) ree 
_ x (6-42 + 16 - 8V2 — 2(2 - v2)V16 — 8V2 — 12 + By2) 


64 + 54 — 36V2 — 32 + 16V2 


_x (22 — 12¥2 - 2(2 - v2)(2)) _x(7-4V2) _ x(7 — 4V2)(43 + 10V2) 
2(43 — 10V2) ~ 43-10V2 (43 — 10V2)(43 + 10V2) 


_ x(301 —172V2 + 70V2— 80) _ x(221 — 102V2) 
7 1849 — 200 Z 1649 


_ x(221 - 102V2) 
a 1649 


T4 = 0.046543490987231234% ae (5) 


Above Eq(5) gives the radius of fourth inscribed circle (see the above figure-6) in terms of edge length x of square 
ABCD. 


Radius of fifth inscribed circle (r,;): Similarly, the radius 7, can be computed by substituting a=, b= 
x/2 and R = x in the above Eq(1) as follows 


abR (ak + DR -ab-2 abR(R — a— b)) 


Cmin = 
me (aR + bR + ab)? — 4abR? 
> 1s 
x(221-—102V2) x _ [ x(221-102y2) x x(221—102V2) x , x(221—102V2) x x(221-102V2) x 
1649 mee 1649 eR ee 1649 oe 1649 a ae 1649 oe 


2 
x(221 — 102V2) x x(221-102V2) x x(221-102V2) x 
( 1649 as ean 1649 A eee 1649 oper 


x (69649 — 450842 + 364429 — 168198V2 — 2(221 — 102V2)V/364429 — 168198V2 — 139298 + 90168V2) 
2719201 + 626841 — 405756V2 — 728858 + 336396V2 


x (434078 — 213282V2 — 2(221 — 102V2)(459 — 85V2)) x(196520 — 82076y2) 
2617184 — 69360V2 ~—-4624(566 — 15v2) 


_ x(170—71V2) _ x(170 — 71V2)(566+15V2) _ x(47045 — 18818y2) 
~ 4(566—15v2) 4(566 — 15V2)(566 + 15V2) — 639812 


_ x(47045 — 18818V2) 
i 639812 


= 0.03193489522432073x% seas (6) 


Above Eq(6) gives the radius of fifth inscribed circle (see the above figure-6) in terms of edge length x of square 
ABCD. 


3.1 Recurrence Relation 


Similarly, the radius rg of sixth inscribed circle can be computed by substituting a = 15, b=x/2 andR =x in 
the above Eq(1) and so on. Every time the case becomes similar to the case of three externally touching circles 
which internally touch the fourth circumscribing circle which is analysed in the above section 1.1. This way we 
can compute the radius of any nth blue circle inscribed in the bounded plane region (see above figure-6). It’s 
worth noticing that the value of radii b = x/2 and R = x are constant (while radius a is variable) in each case 
and the radius 7, of nth inscribed circle depends on the value of radius 7,_, of (n-1)th inscribed circle. 


It can be inferred from the above derivations that the radius 7,,,, of (n+1)th circle inscribed in the bounded region 
can be computed by substituting a = %,, b = x/2 and R = x in the above Eq(1) as follows 


abR (ar + bR—ab —2,/abR(R —a— b)) 
al (aR + bR + ab)? — 4abR? 


2 


2 
x 
(Sem + =) — 2x37, 


5x ta (Ty +x — 2,/™m(x — 2m)) 


x2 
ae (3m, + x)? — 2x3n, 


Ty +X y_ — 2 pV xXT_ — 2772 . 
~ Tay = x| ——— oii 
tet 9r,,2 — 2xr, + x? (7) 


Where, ry = 7 Vx>0,xER, nEN &> yyy =f(ra) 


The above Eq(7) is the governing equation i.c. recurrence relation for 2D packing of circles to analytically 
compute the radius of (n+1)th inscribed circle by substituting the value of radius of nth inscribed circle (i.e. blue 
inscribed circles as shown in the above figure-6). It is worth noticing from above Eq(7) that next term 7,4, is a 
function of the previous term 7, for given side x of the square. 


3.2 MATLAB code-1 


It is usually very difficult to compute the radius of nth (blue) inscribed circle when n is very large therefore using 
the above recurrence relation i.e. Eq(7), a for-loop can easily be used to develop a simple MATLAB code, to 
compute the radii of n number of inscribed circles when the edge length x of square is given, as follows 


clc; 

clear all; 

close all; 

x=input('Enter the edge or side of square (in required unit) = '); 
n=input('Enter the number of circles to be packed, n = '); 


disp('Column matrix of radii of n no. of packed/inscribed circles, r(i)(units)=' ); 
r(1)=x/4; 
disp(r(1)); 
for i=1:1:n-1 
r(1)=x/4; 
*#HCR?s recurrence relation to compute radii of inscribed circles 
r(i+1)=double(x*((r(i) )*2+x*r (i) -2*r (i) *sqrt(x*r (i) -2*(r(i))%2))/(9* (ri) )%2- 
2*x*r(i)+x%*2)); 
A=[r(i+1)]5 
disp(A); 
end 


Above MATLAB code of recurrence relation i.e. Eq(7) requires two input values 1) side x of square and 2) number 
n of circles to be inscribed/packed to compute the radii of all n number of circles packed or inscribed in the plane 
region bounded by given square, semi-circle and quarter circle as shown in the figure-5 above. 


3.3 Results 


After entering the values of side of unit square i.e. x = 1 and number of circles to packed/inscribed i.e. n = 100, 
the following are the radii of all 100 packed/inscribed as tabulated below. 


Inputs 


Enter the edge or side of square (in required unit) = 
1 

Enter the number of circles to be packed, n = 

100 


Outputs 
Column matrix of radii of n no. of packed/inscribed circles, r(i)(units)= 


Table-1: Ratio of radius of packed/inscribed circle to the side of square (n=100 circles) 


@.2500 @.0014 3.7801e-04 1.7120e-04 
Q@.1277 8.0013 3.6374e-04 1.6681e-04 
0.0732 @.0012 3.5025e-04 1.6258e-04 
0.0465 8.0012 3.3751e-04 1.5851e-04 
@.0319 0.0011 3.2544e-04 1.5460e-04 
@.0232 8.0010 3.1401e-04 1.5082e-04 
@.0176 9.4996e-04 3.0318e-04 1.4719e-04 
@.0137 8.9405e-04 2.9289e-04 1.4368e-04 
0.0110 8.4293e-04 2.8312e-04 1.4030e-04 
0.0091 7 .9607e-04 2.7383e-04 1.3703e-04 
@.0076 7.5301e-04 2.6499e-04 1.3388e-04 
0.0064 7.1336e-04 2.5657e-04 1.3083e-04 
0.0055 6.7675e-04 2.4855e-04 1.2789e-04 
0.0048 6.4289e-04 2.4089e-04 1.2505e-04 
0.0042 6.1151e-04 2.3359e-04 1.2230e-04 
0.0037 5 .8236e-04 2.2661e-04 1.1964e-04 
0.0033 5.5526e-04 2.1994e-04 1.1706e-04 
@.0029 5.3000e-04 2.1356e-04 1.1457e-04 
@.0026 5 .0643e-04 2.0745e-04 1.1216e-04 


@.0024 4.8439e-04 2.0161e-04 
8.0022 4.6376e-04 1.9600e-04 
@.8020 4.4442e-04 1.9063e-04 
8.0018 4.2627e-04 1.8547e-04 
0.0017 4.0920e-04 1.8052e-04 
8.0015 3.9314e-04 1.7577e-04 


1.0982e-04 


1.0755e-04 


1.0536e-04 


1.0323e-04 


1.0116e-04 


9.9157e-05 


It’s worth noticing that each value of above Table-1 shows the ratio of radius of a packed/inscribed circle to the 
side of square (i.e. taken as unity). Therefore, the exact value of radius of any of packed/inscribed circle of any 
order (i.e. up to n = 100) can be obtained by directly multiplying given side of square x with the corresponding 
value of ratio as obtained in the above Table-1 (each column has 25 values of ratio) for 100 packed circles. 


Note: The output values of radii of inscribed circles are obtained in form of a column from the above MALAB 
code-1 which have been rearranged in sequence in form of above Table-1 with four columns. 


4. Packing of circles in the plane bounded by two semi circles and a quarter circle 


Consider a square ABCD of each side x which has two semi-circles with 
diameter AB & AD and a quarter circle with centre at the vertex B (see 
figure-7). 


Now, consider the circles with centres C,,C>,C3,...... C,, and radii 
Typ 104 1 3y ve vee T, Yespectively snugly fitted or packed in the plane region 
bounded by two semicircles and a quarter circle inside the square ABCD (as 
shown in the figure-7). Let’s first find the radius 7, of first packed/inscribed 
circle then the radii of all other inscribed circles 7,73, ... ... T, can be easily 
computed by using mathematical relation i.e. Eq(1) similar to the above 


section 3. We have to find out the radius of nth inscribed circle using 


recurrence relation i.e. Eq(7) as derived in the above section 3. 


Now, join the centre C, of first packed/inscribed circle to the centre E & G of 


B 


>| 


Figure-7: Infinite no. of small 
circles can be packed/inscribed in 
the region bounded by _ two 
semicircles and a quarter circle 
inside the square ABCD 


semi-circles and to the vertex B of square ABCD. Drop the perpendiculars C, F and C,H from the centre C, to the 


sides AB and AD respectively (as shown in the figure-8 below). Let C,F = y. 


In right AC, FB, We have, C;F = y and C,B = x —7, (see figure-8 below). Now using Pythagorean Theorem 


as follows 


(C,B)* = (C,F)* + (FB)? 


=> FB=,/(C,B)? — (C.F)? 


FB = V(x—1)?-y 


In right AC, FE, We have, C,F = y, CyE = na T, (see figure-8) 


Now using Pythagorean Theorem as follows 


(C,E)? = (CiF)? + (FE)? > FE= ¥(C,E)? — (C,F)? 


2 


FE = (S+n) — y? 


d 


EF E 
x 


From figure-8, we have B 


“Ke > 

FB=FE+EB © (FB)? =(FE+ EB)? 
Figure-8: A circle of radius 7; and 
5 center C, is the first largest circle 


inscribed in the region bounded by two 
Jan? —y? a x ‘ > .* semi-circles and a quarter circle inside 
> ( (x-1%)?-y ) = (S+n) =) a) the square ABCD. 


x 2 x : ° 
=> (5-3n) = (5+) ~y? ( LHS > 0, RHS > 0) 
2 2 


x 
a ton? — 3xr = 2 Peay 


y? = 4xr, — 8r,? Sirdetentens C41) 
In right AC; HG, We have, 


x x x 2 x 
C.G=5—n, CH= AF = AE~FE=>— (S+n) —y? and HG = HA~GA =(,F~GA=y—> 


Now using Pythagorean Theorem (see figure-8 above) as follows 


ll 


(C,G)* = (C,H)? + (HG)? 


2 
2 x x 2 xX\2 


2 4p 8) |G) = ae 7) 


x? ae ia 
—+r?%—-xn =—+—47r724+x4- +n) —y?+y? tena 
4 4 4 

= + 2r, y= | =+r) (- x #0) 


= (5+2n-y) = (F+n) ~y? (: LHS > 0, RHS > 0) 


2 2 


x x 
qt an’ ty? + 2xn — Any — xy = atm txn-y? 


2y? —(x+4n)y + xr, + 3r,? =0 


2(4xr, — 87,7) — (x + 4rn)y +xr, +37? =0 (Setting value of y* from Eq(4.1)) 
9xr, — 137,? 
= . (4.2 
x+4r, (42) 
Now, substituting the value of y from Eq(4.2) into above Eq(4.1) as follows 
9xr, — 137,2\* 
4) = 4xr, — 81,2 
x+4n 
ry? (9u— 135)? Ana 2h G+ 4h) 
(x + 4r,)2 ~ 
r{r, (81x? + 1697,7 — 234xr,) — (4x — 87,)(x? + 167r,7 + 8xr,)} = 0 (- x+4r, #0) 
81x2r, + 1697,3 — 234x7,? — 4x3 — 64xr,? — 32x77, + 8x77, + 1287,3 + 64x7,? = 0 (- 7, #0) 
2977r,3 — 234x172 + 57x27, — 4x3 = 0 
13 ny? TY _ er 3 
297 ) ~ 234 (<) +57 (=) —~4=0  (Dividingbyx?) __.....(4.3) 


The above Eq(4.3) is a cubic equation in terms of which can be easily solved by numerical techniques which 


give the following solutions 


4 
"33 


r, +11 
x 3'3 


It is to be noted that the region bounded by two semicircles and a quarter circle inside the square ABCD (see the 
above figure-7) is less than the region bounded by semicircle, quarter circle and the square ABCD (see the above 
figure-4). This implies that the radius 7, of first inscribed circle in this case (section 4) must be less than the radius 


1 = - obtained in the above section 3. Therefore, We select the real roots of above cubic equation (4.3) such that 


= <jZorn< - . After applying this condition, we obtain a unique solution as follows 


Tr, 4 4x x x 
x 33 og ee) 
_ 4x (8) 
“T1= 33 a 


The above Eq(8) gives the radius of first largest packed/inscribed circle (as shown in the above figure-7) in terms 
of edge length x of square ABCD. 


Coordinates of the centre (C,): The coordinates of the centre C, of first inscribed circle w.r.t. the vertex A 
(assumed to be the origin) of the square ABCD (see above figure-8) can be determined as follows 


4x 4x2 
di =CF= _ 9xr, 1372 9x" 3q— 13 (g3) _ 980x 20x 
y — coordinate = C,|F =y= = _ an ee Pea 


x x 2 x x 4x\? 20x,” 
x — coordinate = C,H = AF = AE — FE = > — (S+n) =ye Soo ( + ) -( ) 


2 2 33 33 
= 4x 
4 
Therefore, the coordinates of the centre C, of first inscribed circle w.r.t. the vertex A are (=, =) (see the figure- 


8 above). Similarly, the coordinates of other inscribed circles can be computed after computing their radii. 


Radius of second inscribed circle (12): Now, draw the complete circles by extending the semi-circle with centre 
E and quarter-circle with centre B (as shown in the figure-9 below). It is worth noticing that three externally 
touching circles with centres C, ,C, and E touch the circumscribing circle (i.e. extended form of quarter circle in 
the above figure-7) with centre B. This case is similar to the case of three externally touching circles which touch 
the forth circumscribing circle analysed in the above section 1.1 (see the above figure-1). 


Figure-9: The semi-circle and the quarter circle are extended to draw 
complete circles with centers E & B and radii x/2 & x respectively. 


Now, the radius 7, of the second inscribed circle with centre C, touching the circles with centres C,, E & B and 
the radii 7, = 4x/33, x/2 & x respectively, is found out by substituting a = 17, = 4x/33, b=x/2 andR= 
x in the above Eq(1) (i.e. minimum value Cin) as follows 


abR (aR + DR - ab—2 abR(R — a— b)) 


nS (aR + bR + ab)? — 4abR? 
4x x 4x x 4x x 4x x 4x x 
332 «(33 x+7°X— 33°97 — 2/3379 ( 33 >) 
> n= 


2x3 (ae a) 68x° 


2, BINGE: - B82} e  BEe de 
~ 745 y2\2 gx4 969x457 
( 22 ) ~ 33 667 
4x 


T2= 57 0.07017543859649122 x 


The above Eq(9) gives the radius of second inscribed circle (as shown in the above figure-7) in terms of edge 
length x of square ABCD. 


Radius of third inscribed circle (r3): Similarly, consider three externally touching circles with centres 
Cz ,C3 and E which internally touch the largest circle with centre B (as shown in the figure-7 above). This case is 
again similar to the case analysed in the above section 1.1. Thus, the radius 73 of circle with centre C3 can be 
computed by substituting a = 1r,, b = x/2 and R = x in the above Eq(1) as follows 


abR (ar + bR — ab — 2,fabR(R —a— b)) 
min = (aR + bR + ab)? — 4abR? 


4x xX. Ax 4X, _4x xX _ 5 Ax XxX. (x- 2-3) 
57.2 NSE oe eS BF 2 ke G7 2 


eae 4 dx oe. 4 
x. x. X XV yg FX XL? 
ee ey ee 
2x3 (61x? 14x? BGxe 
_SO7\TI4 57 )_ tat _ 4 
~ 169x2\2 Bxt 2937x489 
(Sir ~ 37 114* 
4x 
T3 = By © 0.0449438202247191 x oases (10) 


The above Eq(10) gives the radius of third inscribed circle (as shown in the above figure-7) in terms of edge length 
x of square ABCD. 


Radius of forth inscribed circle (74): Similarly, the radius 7, can be computed by substituting a = 73, b = 
x/2 and R = x in the above Eq(1) as follows 


abR (ar + bR — ab —2,/abR(R —a— b)) 
ee (aR + bR + ab)? — 4abR? 


eS 4 4 2 4 
x x xX XV 4. FX XV 2 
(gap eta g) Ago oe 
2x3 Ont’ Taxes 228x° 
eo i178 89) Ger _ 4x 
~ 7101x2\2 gx4 7353x4129 
( 178 ) ~ Bo 178? 
4x 
14 = G59 © 0-031007751937984496x oe (11) 


The above Eq(11) gives the radius of fourth inscribed circle (as shown in the above figure-7) in terms of edge 
length x of square ABCD. 


Radius of fifth inscribed circle (7): Similarly, the radius 7, can be computed by substituting a=, b= 
x/2 and R = x in the above Eq(1) as follows 


abR (ar +bR—ab —2,/abR(R —a— b)) 
a (aR + bR + ab)? — 4abR? 


Ax xX. Ax | 4%, _ AX xX _ 5 Ax XxX. (<-t5-3) 
129° 20" jog. Ok 720 129: 2° “99 2 


> Ts = 4 4 2 
xX x xX XxX xX XxX 
Gua o eae) a 
De (Bake 22K" 356x° 
129258 129) = gsgxr 4x 
~444y%2\2 gx4 —«15753x* 177 
( 258 ) ~ 729 258? 
4x 
rs = ——= © 0.022598870056497175xX antec (12) 


177 


The above Eq(12) gives the radius of fifth inscribed circle (as shown in the above figure-7) in terms of edge length 
x of square ABCD. 


Similarly, the radii 7g, 77, 7g, .....% of all other inscribed circles can be found out. 


It is worth noticing that the value of radius of each inscribed circle comes out in form of p/q whose numerator is 
always constant 4x (x is constant side of square) but their denominators; 33, 57, 89, 129, 177, ........ form a 
series in which the differences between consecutive terms form an Arithmetic Progression (A.P.). Let this 
unknown series of the denominators be expressed in form of the sum of n number of terms as follows 


Se 93 G7 89 129-177 es aly Es (1) (T, is nth term) 
Sy = 33 4+574+ 8941294177 4 wee $Ty (2) (Rewriting the series) 


Now, subtracting the above series (2) from (1) (i.e. taking difference of corresponding terms of above two series), 
we obtain 


0 = 334+ 244 32 +40 +484. cece ees In (Total (n + 1) terms on RHS) 
Ty — 33 = 244 32 + 40 + 484. eee cee ee ee cee cee ee (Total (n — 1) terms on RHS) 


The series on the RHS is an Arithmetic Progression (A.P.) with first term a = 24, common difference d = 8 and 
total (n — 1) terms. Now, taking the sum of all (n-1) terms of above A.P. (on RHS) as follows 


(n — 1) 
2 


E332 (2:24+(n-1-1):8) =4n-DM+4 


T, = 33+4n? +12n—16 
T, = 4n* +12n+17 


Since the numerator of p/q form of radius of inscribed circle remains constant equal to 4x but the denominator of 
nth p/q form i.e. radius r, of nth inscribed circle is T, = 4n?+12n+17 . Therefore, following the 
symmetry/pattern in the denominators of p/q form of radii 11, 72,73,1%4,15 .......% Of inscribed circles (as derived 
in above section 4), the radius 7, of nth inscribed circle is given as follows 


4x 4x 
"TT, 4n2+12n+17 


4x 


= —>—_—_——\{" Si tetateer UZ 
4n2 +12n+17 oe 


as 


The above Eq(12) is the generalized formula to directly and analytically compute the radius 7, of nth circle 
packed/inscribed in the plane region bounded by two semicircles and a quarter circle inside the square of given 
edge length x (as shown in the above figure-7). 


Now, substituting the different values of natural number n = 1, 2,3,4...... in the above generalized formula i.e. 
Eq(12), the radii of inscribed circles are obtained as follows 


4x 4x 4x 4x 5 4x a 4x - 4x 2 4x 
gg. a7? gg" 190" ayy? 2as" 907 = gear 


NY — 


The recurrence relation i.e. Eq(7) is equally applicable in this case to find the radius of any inscribed circle as 
follows 


aD) 


9r,2 — 2xr, + x? 


Tye +X _ — 2 pV XT _ — 27 n?2 
Tau = x| ————____—_ 


Where, r= = Vx>0,xER, nEN & ras =f(rn) 


However a simple MATLAB code has been developed using recurrence relation i.e. Eq(7) to compute the radius 
of all n number of circles packed/inscribed in the plane region bounded by two semicircles and a quarter circle 
inside a square of given edge length, as follows 


4.1 MATLAB code-2 


Using the above recurrence relation i.e. Eq(7), a for-loop can easily be used to develop a simple MATLAB code 


clc; 

clear all; 

close all; 

x=input('Enter the edge or side of square (in required unit) = '); 
n=input('Enter the number of circles to be packed, n = '); 


disp('Column matrix of radii of n no. of packed/inscribed circles, r(i)(units)=' ); 
r(1)=(4*x)/33; 
disp(r(1)); 
for i=1:1:n-1 

r(1)=(4*x)/33; 
*#HCR?s recurrence relation to compute radii of inscribed circles 
r(it+1)=double(x*((r(i))*24+x*r(i)-2*r(i)*sqrt(x*r(i)-2*(r(i))%2))/(9* (ri) )%2- 
2*x*r(i)+x%*2)); 
A=[r(it1)]; 
disp(A); 
end 
Above MATLAB code of recurrence relation i.e. Eq(7) requires two input values 1) side x of square and 2) number 
n of circles to be inscribed/packed to compute the radii of all n number of circles packed or inscribed in the plane 


region bounded by two semi-circles and a quarter circle inside the square as shown in the figure-7 above. 


4.2 Results 


After entering the values of side of unit square i.e. x = 1 and number of circles to packed/inscribed i.e. n = 100, 
the following are the radii of all 100 packed/inscribed as tabulated below. 


Inputs 


Enter the edge or side of square (in required unit) = 
1 

Enter the number of circles to be packed, n = 

100 


Outputs 
Column matrix of radii of n no. of packed/inscribed circles, r(i)(units)= 


Table-2: Ratio of radius of packed/inscribed circle to the side of square (n=100 circles) 


@.1212 0.0013 3.6255e-04 1.6644e-04 
8.0702 @.0012 3.4913e-04 1.6223e-04 
0.0449 0.0011 3.3645e-04 1.5817e-04 
8.0310 0.0011 3.2444e-04 1.5427e-04 
@.0226 8.0010 3.1306e-04 1.5051e-04 
@.0172 9.4496e-04 3.0227e-04 1.4688e-04 


0.0135 8.8948e-04 2.9203e-04 1.4338e-04 
0.0108 8.3875e-04 2.8231e-04 1.4001e-04 
0.0089 7.9224e-04 2.7306e-04 1.3676e-04 
0.0074 7.4948e-04 2.6425e-04 1.3361e-04 
0.0063 7.1010e-04 2.5587e-04 1.3058e-04 
0.0054 6.7374e-04 2.4788e-04 1.2764e-04 
0.0047 6.4010e-04 2.4025e-04 1.2481e-04 
0.0041 6.0892e-04 2.3298e-04 1.2207e-04 
0.0036 5.7996e-04 2.2603e-04 1.1941e-04 
@.0032 5.5302e-04 2.1938e-04 1.1685e-04 
@.0029 5.2791e-04 2.1303e-04 1.1436e-04 
@.0026 5.0448e-04 2.0694e-04 1.1195e-04 
@.8024 4.8257e-04 2.0112e-04 1.0962e-04 
@.0022 4.6205e-04 1.9553e-04 1.0736e-04 
8.0020 4.4282e-04 1.9018e-04 1.0517e-04 
8.0018 4.2476e-04 1.8504e-04 1.0305e-04 
@.0017 4.0779e-04 1.8011e-04 1.0099e-04 
0.0015 3.9181e-04 1.7537e-04 9.8988e-05 
8.0014 3.7675e-04 1.7082e-04 9.7047e-05 


It’s worth noticing that each value of above Table-2 shows the ratio of radius of a packed/inscribed circle to the 
side of square (i.e. taken as unity). Therefore, the exact value of radius of any of packed/inscribed circle of any 
order (i.e. up ton = 100) can be obtained by directly multiplying given side of square x with the corresponding 
value of ratio as obtained in the above Table-2 (each column has 25 values of ratio) for 100 packed circles. 


Note: The output values of radii of inscribed circles are obtained in form of a column from the above MALAB 
code-2 which have been rearranged in sequence in form of above Table-2 with four columns. All the values (i.e. 
obtained by recurrence relation) of Table-2 can be easily verified and validated by comparing them with the values 
obtained by above analytical formula i.e. Eq(12). 


5. Packing of circles inside the sector of given radius and aperture angle 


Consider a sector OAB consisting of two radii OA and OB each of length R 
and an intercepted arc AB such that the central angle exerted by arc AB at the 
centre O is 6, (as shown in the figure-10). 


Now, consider n number of circles with centres C,,C3,C3,...... C,, and radii 
Tis Voy Fay aie T, respectively packed/inscribed in the sector OAB (as shown in 


the figure-10). Let’s first find the radius 7, of first packed/inscribed circle then 
the radii of all other inscribed circles 73,73, ...... 7 can be easily computed by / ) 


using recurrence relation. e———_— Rk 


: . . . Figure-10: Infinite no. of small 
Join the centres C; and C, and to the centre O of sector OAB. Drop the Gircles with centers C4, Cy, Cy 


perpendiculars C,M and C,N from the centres C, and C2 to the radii OA and can be packed/inscribed in the 


OB respectively (as shown in the figure-11 below). secon Ooms 
In right AC, MO, We have, OC; = R—1,, C,M = 7, (see figure-11 below) 
‘ C,M _ 94 "% 
sin ZC,0M = 0c, > sin a 7 
ee Oy, 1+sin 2) = Resin? 
sin> —nsin> = 7 > r,( sin - sin > 
R sin 
sem (V0<0,<7) Bicneawi( 13) 

1+ sin 7 


The above Eq(13) gives the radius of first inscribed circle (as 
shown in the above figure-10) in terms of radius R and central 
angle 6, of sector OAB. 


In right AC,NO, We have, C,N = 7, OC, = R —79. 


Let ZBOD = 6) (see figure-11). 


in ZC.ON 2N . 9, TY) 
> sin Z =— > sin—= 
‘ Oc, 2 R-1% 
"% 
one " R x 
sin ria 7 4] 
R( ad 3) I= R * Figure-11: Infinite no. of small circles with 
centers C,,C2,C3... can be packed/inscribed in 
1% . : . ; the sector OAB. 
Where, x = ; is Some variable depending on the value of radius 
r2 for given radius R of sector AOB. 
02 1 , > 92 ( 0<86 9 0 6, 9 1 02 0) 
=> cos—= — sin* — ~0<0.<0,<01 SO 0<><><=>5 cos—> 
2 2 aii a 2 


In AC, OC,, We have, OC; = R —17,, OC, = R—71, and C,C, = 7 +7, . (see above figure-11). 


A 


2 2 2 


ZAOB ZBOD 6, 6, 
=> ZC,0C, = ZC,0B — Z2C,0B = = 


Now using cosine rule in AC, OC, as follows 


(OC,)? + (OC,)? — (C,C,)? 


ZC,0C, = 
Page 2(0C,)(OC,) 


0, 8 (R-1n)*+(R-%)*?-(% +7)? 
cos ( - ) = 


2 2 2(R —17,)(R — 1) 
0, On 8, 6, R*-Rr,— Rr, - 7 
cos — cos sin sin = ER =A) 


ign ke 6 . 0 3 : 
Snow, substituting the values of cos = , sin = and 7, in the above equation as follows 


Rsin@ Rsin@ 
R? — -Rr - Tr. 
6, V1—2x f) 1+sin+ 1+sin+ 
cos . 1 +sin> 1 = 4 7 4 
mi - Rsin=> 
= 6 (R—-1) 
1+sin> 
6, o5Gia. Hoe oe g, (1 + sin — sin — 3 - BsinZ — Bsin) 
1—2xcos> + x sina 1+ sin= 
1—-x R2 Oe- -.aBy T 
——z (1+ sin - sin) (1-F) 
1+ sin= 
Tonos oe een (1 —x—2x sin) 
Z ges 2 (. 2 =x) 
1-x 1-x R 
0, . 0, . 0, 
POPES ge ese Le ee ee 
1-x 
0, 0, 
T= x cos—+(1+3sin)x-1=0 (vx-140>241Vv 7, #R) 


pr (1435 2) 
x COS > = sin>}* 
Taking the squares on both the (positive) sides of the above equation as follows 


6,\° 6, 
( T= 2x cos) =(1-(143sin3)x) 


24 k 0,\7 2 i 0; 
(1 — 2x) cos >= 1+ 3 sin x°-—2 1+ 3sin> x+1 


2 


64\? 6g 6 6 
(1 +3sin2) x?—2(1+3sin cos’) x +1- cos? = 0 
2 


0 0 0 0 
(1 + 3 sin) x? — 2sin (3 + sin) x + sin? = =0 


Solving the above quadratic equation for x using quadratic formula as follows 


_ (-2 sin (3 + sin) + (-2 sin (3 + sin %)) —4 (1 + 3 sin Sy sin2 a 


2 (1 #3 sin) 


2 
. 0, . 0, . 0, ie . 2 94 
ry _ 2sin (3 + sin) + 2sin 8 — 8sin 2 G _") 
Rr 6.\2 ney: 
2 (1+ 3sin) 
. 0, . 0, . 0, 0, 
_ Rsinp (3 + sin Fp) + 2V2sincos-F 
TSS rR 
Oy 
(1+ 3sin3) 
R (sin? +3 sin + V2sind,) 
i) a 


(1 +3 sin) 


Case-1: Taking positive sign 


R (sin? a +3 sin + V2sin ,) 
n= Sn fae 
(1 +3sin 3) 
The values of radii obtained from the above expression of 7, when used as a recurrence relation have been 
compared with the values of radii obtained from the drawings and measurements. It has been verified that the 
above expression of the radius ry gives erroneous values of radii of the packed/inscribed circles (i.e. n number of 
circles) in the sector because the recurrence relation obtained from the above value of Tr, gives the increasing 
values of radii of inscribed circles with decreasing size which is not possible. Therefore, this case is discarded. 


Case-2: Taking negative sign 


R (sin? a +3 sin —v2sin ,) 
ie 


ei Ne 
(1 +3sin 4) 
It has been verified that the above expression of 72 when used as a recurrence relation gives the decreasing values 
of radii of inscribed circles (i.e. n number of circles) which are exactly same as obtained from the drawings and 
measurements. Hence this case is valid and the above value of 72 is acceptable. 


R (sin? a +3 sin — 2 sin 61) 
a eee seated: ke Ay 


(4 +3sin By. 


The above Eq(14) gives the value of radius r, of second circle inscribed in the sector OAB (as shown in the figure- 
11 above). 


Now, substituting the value of 7, from Eq(14) in the expression of sin 2 (i.e. derived above) to obtain the value 


angle 8 as follows 
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R (sin? A143 sin —y2sin 6;) 


2 
2 
0, 1 (1 + 3sin &) sin? a +3 sin — v2sin6, 
sin = a = 

2 2G R (sin? 2 + 3 sin — y2sind,) 1+ gsin? 2 + 3sin% + V2sin 8, 

= ne Sa. re 

(1 +3sin 3) 

sin? 21 + 3sin 91 — y2sin 0, 

0, =2sin1 be cideseatt lo) 


V2 sin, + gsin? 91 + 3 sin 32 + 1 
The above Eq(15) gives the value of angle 8, exerted by the second inscribed circle at the centre O of the sector 
OAB (as shown in the figure-11 above). 


Now, consider the new sector BOD of radius R and 
central angle 0, which inscribes the second circle of 
radius 7, and centre C, and the third circle of radius 
Tr, and centre C, (as shown in the figure-12). This 
case becomes similar to the precious case in the 
above figure-11 in which the sector OAB of radius R 
and central angle 0, inscribes the circles with centres 
C, and Cy. 


Therefore using the symmetry in the sectors OAB 
and OBD (see figure-11 & figure-12), the radius 73 z 
of third packed/inscribed circle with centre C3 can be Oo R D 


computed by substituting 6, = 6, in the above Figure-12: The sectors BOD and BOE are similar to the 


Eq(14) as follows original sector OAB of radius R and central angle 0,. 


R (sin? 2 + 3 sin 2 — V2 sin 6, ) 


(1 + 3 sin 2)" 


13 = 
The angle 03 subtended by the third packed/inscribed circle with centre C3 at the centre O of the sector OBD (i.e. 
sector OAB ) (see figure-12) can be computed by substituting 8, = 02 in the above Eq(15) as follows 


7) 7) 
in? + in 2 — V2 sj 
‘ sin 7) 3 sin 7) 2 sin 02 


63 =2sin— 6 7 
V2 sin 0, + 8sin? + 3sin= +1 

Similarly, consider the sector OBE of radius R and central angle 63 which inscribes the third circle of radius r3 
and centre C3. The radius 7, and the subtended angle 0, of the fourth packed/inscribed circle are obtained by 
substituting 8, = 83 in the above Eq(14) and Eq(15) as follows 


R (sin? 3+ 3 sin 9 — V2 sin 05) 


2 
i= a2 
in 
(1+ 3 sin 3) 
sin? 2 + 3 sin — y2sin 6, 
V2 sin 63 + 8sin? $+ 3sin=+1 
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Similarly, the radii 75,176,177 «++ and exerted angles 05,606,607 ...... of all other circles packed/inscribed in the 
original sector OAB of radius R and the cetral angle 0, can be determined by substituting 0, by the previous value 
in the above Eq(14) and Eq(15). 


It can be concluded that the radius 7,,, and subtended angle 6,,, of say (n + 1)th packed/inscribed circle 
depends on the previous value of angle 6, i.e. angle exerted by nth packed/inscribed circle. Therefore, the radius 
Tp+1 and subtended angle 6,,,; of (n + 1)th packed/inscribed circle are determined by substituting 6, = 6, in the 
above Eq(14) and Eq(15) as follows 


R (sin? Bn +3 sin — V2 sin 8) 
Tn+1 = 2 ous sen cee eee (16) 


(1+3sin 2) 


sin? Sn + 3 sin 5 — V2sin 8, 


sin 6, + 8sin* —- + 3 sinz=- + 
V2 sin 6, + 8si 2 on 3 si on 1 


1 


On41 = 2sin7 in pecan (LZ) 


O44 
a v 0<4<T7, n€N = Tnii = f(On), 9n+1 = g(On) 


The above recurrence relations i.e. Eq(16) and Eq(17) can be used to determine the radius and the angle subtended 
by any of infinite number of circles packed/inscribed in the sector of radius R and central angle 6, (refer to the 
above figure-10). 


It is interesting to note that the radius and the subtended angle (i.e. obtained from Eq(16) & Eq(17)) of 
packed/inscribed circle decrease dramatically i.e. their decay is much faster than the exponential decay. In other 
words, the values of radius (7;) and the subtended angle (9;) approach to zero faster while packing a few number 
of circles in the sector of finite radius and central angle. 


A simple MATLAB code has been developed using the above recurrence relations i.e. Eq(16) and Eq(17) to 
compute the radii and the subtended angles of all n number of circles packed/inscribed in the sector of given radius 
and central angle, as follows 


5.1 MATLAB code-3 


Using the above recurrence relations i.e. Eq(16) and Eq(17), a for-loop can easily be used to develop a simple 
MATLAB code 


clc; 

clear all; 

close all; 

R=input('Enter the radius of original sector, R (in required unit) = '); 
t1i=input('Enter the central angle of sector, theta_1 (in degree) = '); 
n=input('Enter the number of circles to be packed, n = '); 
r1=double((R*sind(t1/2) )/(1+sind(t1/2))); 

A=[r1 t1]; 


disp('(nx2) matrix of radii & angles subtended by inscribed circles, [r(units) 
theta(degrees)] ='); 
disp(A); 
for i=1:n-1 
t(1)=(t1/2)*(pi/18@) ; 
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t(it+1)=asin(((sin(t(i)))*2+3*sin(t(i))- 
sqrt(2)*sin(2*t(i)))/(8*(sin(t(i)))*2+3*sin(t(i))+sqrt(2)*sin(2*t(i))+1)); 
r(it+1)=double((R*((sin(t(i)))*2+3*sin(t(i))- 
sqrt(2)*sin(2*t(i))))/(1+3*sin(t(i)))%2); 


B=[r(i+1) (360/pi)*t(i+1) ]; 


disp(B); 


end 


5.2 Results 


Using the above MATLAB code, the values of radii 7; (in given unit) and the corresponding values of angles 6; 


(in degree) subtended by 15 circles packed/inscribed in the sectors of unit radius (R = 1) and central angles 6, = 
30°, 45°, 60°, 90°, 120°, 135°, 150°, are obtained and arranged in the Table-3 below. 


Table-3: Radii and angles subtended by the circles packed/inscribed in a sector of unit radius. 


i N% 0; N% 0; N% 0; r% 0; rj 6; % 6; T% 6; 

1 0.2056 30.0000 0.2768 45.0000 0.3333 60.0000 0.4142 90.0000 0.4641 120.0000 0.4802 135.0000 0.4913 150.0000 

2 0.0432 3.1759 0.0638 7.8185 0.0840 10.5288 0.1239 16.2602 0.1640 22.6285 0.1845 26.1597 0.2056 30.0000 

a 0.0077 0.8882 0.0116 1.3420 0.0155 1.8077 0.0238 2.7944 0.0329 3.8947 0.0378 4.5073 0.0432 5.1759 

4 0.0013 0.1524 0.0020 0.2302 0.0027 0.3102 0.0042 0.4795 0.0058 0.6683 0.0067 0.7734 0.0077 0.8882 

5 0.0002 0.0261 0.0003 0.0395 0.0005 0.0532 0.0007 0.0823 0.0010 0.1147 0.0012 0.1327 0.0013 0.1524 

6 39.08 x 10-° 0.0045 59.14 x 107° 0.0068 79.65 x 107° 0.0091 12.32 x 10-5 0.0141 L7A7x L0r* 0.0197 19.86 x 10-5 0.0228 22.81 x 1075 0.0261 

7 67 410°" 76.97 x 10-5 10.18 x 10-¢ 0.0012 13.62 x 107° 0.0016 21.11 x 107° 0.0024. 29.49 x 1078 0.0034 34.14 x 107° 0.0039 39.08 x 107° 0.0045 

8 12 x 107’ 13,21 x 10-5 17.%.107 19.95 x 1075 23% 107" 26.88 x 1075 36x 1077 41.55 x 1075 S?%.10/% 57.91 x 10-5 58 x 1077 67.02 x 1075 67x 1077 76.97 x 1075 
9 20 x 10-8 22.66 x 107° 30 x 10-8 34.23 x 10-6 40 x 1078 46.11 x 107° 62 x 1078 71.28 x 10-6 87 x 10-8 99.36 x 107° 10x 10-7 11.5 x 1075 12x 107’ 13.21 x 10-5 
10 34x 107% 38.87 x 10°” 51% 107? 58.73 x 1077 69 x 107° 79.12 %10°7 11x 10-8 12,23'% 10°° 15 x 10-8 17.05 x 10-6 17 x 10-8 19.73 x 10-¢ 20 x 10-8 22.66 x 107° 
iat 58 x 10710 66.7 x 10-8 9x 107° 10.08 x 1077 12 x 107° 13.57 ¥ 10-7 18 x 1079 20.98 x 107” 26 x 1079 29.25 x 1077 30x 107° 33.85 x 1077 34x 107° 38.87 x 107” 
12 10 x 10-10 11.44 x 10-8 is 107% 17.29 x 10-8 20x 10710 23.29 x 10-8 31% 107 36 x 10-8 44x 10710 50.18 x 10-8 51x 10-10 58.08 x 10-8 58x 10710 66.7 x 10-8 

13 17x 1071! 19.63 x 107° 26x 107" 29.66 x 107° 35x 1071 39.96 x 107° 54x 1071 61.77 x 107° 75% 1072 86.1 x 10-° 87x 107" 99.65 x 107° 10% 10-79 11.44 x 1078 
14 29 x 107? 33.69 x 10719 44x 10° 50.9 x 10719 60x 10-4 68.56 x 10719 9%, 107% 10.6 x 10° 13'* 207% 14.77 x 107° 15 X 10-7 17.1x 10° 17% 20-2 19.63 x 107° 
15 50 x 10-8 57.80 x 10-14 76x 10-8 87:32 x10" 10x 10" 11.76 x 10-7° 16x 10°? 18.18 x 10-'° 22'« 1077? 25.35\K07 26x 107 29.33 x 10719 29 10-" 33.69 x 1071° 
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6. Packing of circles in the plane bounded by two identical circles & their common tangent 


Consider two externally touching identical circles with centres A & B, and each of radius R and their common 
tangent PQ. (as shown in the figure-13). 


Figure-13: n no. of circles are inscribed in the plane region bounded by two externally 
touching identical circles with centers A & B and each of radius R, and their common 
tangent PQ. The common tangent PQ can be assumed to be a circular arc of infinite radius. 


Now, consider n number of circles with centres Cy, C2, C3,...... C, and radii 14,172,713)... % Tespectively 
packed/inscribed in the plane region bounded by two externally touching identical circles and their common 
tangent (as shown in the figure-13 above). Let’s first find the radius r,of first packed/inscribed circle then the 
radii of all other inscribed circles 73,73, ... ... 1, can be easily computed by using generalized formula. 


We know that the radius (1) of inscribed circle which externally touches three 
given circles with radii say a, b & c touching each other externally (as shown 
in figure-14) is given by the following generalized formula [1] 


abc 
1 ie 
2Jjabc(a+b+c)+(ab+bc+ca) Ko= 


It’s worth noticing that the common tangent PQ can be assumed to be the arc 
of a circle of infinite radius which externally touches two externally touching 
equal circles with centres A and B (see the above figure-13). Thus this case 
(figure-13) becomes similar to the case (figure-14) of three externally 
touching circles with radiia = R,b = Randc > ©. 


Now, the first circle with radius 7, and centre C, is packed/inscribed in the 

: : ee : : Figure-14: Three externally touching 
plane region bounded two equal circles and an infinite circle (i.e. common circles with radii a,b, and c inscribe a 
tangent PQ). The first packed/inscribed circle with centre C, and radius 1, _ blue circle with radius r. 
externally touches three externally touching circles with radii a = R,b=R 


and c — © (see the above figure-13). 


Therefore, the radius 7, of first packed/inscribed circle (see figure-13 above) is determined by substituting a = R, 
b = R in the above Eq(18) and taking its limit as c > © as follows 


; R«Rec : Rc 
r, = lim = lim 


com 2/R-R-c(R+R+c)+(R-R+R-c+c-R)  Re(2 (E+E 1)+B+2)) 
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R _ R _R 
R 2/(0+0+1)+0+2 4 


Now, consider the second packed/inscribed circle with radius r, and centre C, which externally touches the three 
externally touching circles with centres A,B and C, and radii R, R and r, respectively (see the above figure-13). 


Therefore, the radius rz is determined by substituting a = R, b= Randc = = = in the above Eq(18) as follows 


3 
RRS x R 
TB ~ SR? 3R2 12 
R R RR 
pt us . Aenean + 
2|R-R-G(R+R+G)+(R RR He R) 7) 7) 
> T2= 55 


12 


Similarly, consider the third packed/inscribed circle with radius r; and centre C3 which externally touches the 
three externally touching circles with centres A, B and C, and radii R, R and ry respectively (see the above figure- 


13). Therefore, the radius rz is determined by substituting a= R, b = R andc =17, = = in the above Eq(18) as 


follows 
R R3 
r= ne = 12 ee: 
R R R R SR2 7R2 24 
ae Re ida cLiairar ee 
2] R (R+R+G)t(R-R+R +R) Sts 
R 
> =— 
"334 


Similarly, the radius ™ is determined by substituting a= R, b=R and c=73 = = in the above Eq(18) as 


follows 
R R3 
— RR og = A _R 
R R R R 7R2  13R2 40 
pie soe ; dN ah Fe + 
2 |r R sq(R+R +53 +(R R+R agt5q R) itp 
R 
> =_- 
"4 ~ 40 


Similarly, the radius 7, is determined by substituting a= R, b=R andc=%= = in the above Eq(18) as 


follows 
R R3 
RS ae — 40 -— 
R R RR 9R? , 21R* 60 
2IR-Reao(R+R+ay +(R-R+R wtae R) x0 + 30 
R 
> = — 
"s = 60 


Similarly, the radius rg is determined by substituting a= R, b=R andc=r7 = = in the above Eq(18) as 


follows 
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R R? 
T= == 00 = TRGSETEL 7 ai 
afar B(R+r+f)+(a-n+R-A+h-R) 39 +39 
R 
To=o4 
Similarly, all other radii 17, 7g, 79, .....% can be determined using the generalized formula i.e. the above Eq(18). 


It is worth noticing that the value of radius of each inscribed circle comes out in form of ; whose numerator is 


always constant R (i.e. given radius of equal circles) but their denominators; 4, 12, 24, 40, 60, 84, ........ form a 
series in which the differences between consecutive terms form an Arithmetic Progression (A.P.). Let this 
unknown series of the denominators be expressed in form of the sum of n number of terms as follows 


Sy =44+124 244404 604 84 we Th (1) (T, is nth term) 
Sy = 44+12+244+40+60+4 84.0.0... 47, (2) (Rewriting the series) 


Now, subtracting the above series (2) from (1) (i.e. taking difference of corresponding terms of above two series), 
we obtain 


0=44+84124 164 204 244 20. ee Ty (Total (n + 1) terms on RHS) 
Ty, = 4484124 164 20 4+ 244. oe ee (Total n terms on RHS) 


The series on the RHS is an Arithmetic Progression (A.P.) with first term a = 4, common difference d = 4 and 
total n terms. Now, taking the sum of all n no. of terms of above A.P. (on RHS) as follows 


n 
Tr =72:4+(—-1)-4) 
T, = 2n(n+ 1) 


Since the numerator of form of radius of inscribed circle remains constant equal to R but the denominator of 


nth 7 form i.e. radius %, of nth inscribed circle is T,, = 2n(n + 1) . Therefore, following the symmetry/pattern in 


the denominators of form of radii 11,172,173, 1%4,T5 +. +++T%, the radius 7, of nth packed/inscribed circle is given 
as follows 
ome R R 
1 T, 2n(n+1) 
R 


The above Eq(19) is the generalized formula to directly and analytically compute the radius 7, of nth circle 
packed/inscribed in the plane region bounded by two externally touching identical circles each of radius R and 
their common tangent (as shown in the above figure-13). 


Now, substituting the different values of natural number n = 1,2,3,4...... , in the above generalized formula i.e. 
Eq(19), the radii of inscribed circles are obtained as follows 


R R R R R R RoR 
Ha 79) 8 = 54s ag” Eno Ba. tie tae 


y= 
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7. Closest 2D packing of identical circles in a regular hexagon 


The closest, densest or hexagonal packing of circles over a plane is the 
closest/densest possible 2D arrangement of circles without any overlap such that 
each circle touches exactly six circles in its closest neighbourhood. (as shown 
in the figure-15). We are to derive the generalized formula to compute total 
number of identical circles most densely packed in a regular hexagon, the radius 
of smallest circle and the minimum length of string which encloses the identical 
circles of a given (finite) radius which are the most densely packed in a regular 
hexagon of given side. 


Let’s first define the following known parameters Figure-15: The circle with center O 
and radius r is tangent to 6 identical 
circles the most densely packed in a 


n = number of identical circles tangent to any of six sides of a regular hexagon 
regular hexagon. 


in which the circles are the most densely packed 

r = radius of each of identical circles the most densely packed in the regular hexagon 

The following important unknown parameters are to be determined in terms of known parameters (7,7) 

N = total number of identical circles which are the most densely packed in the regular hexagon 

N, = total number of voids (i.e. bounded regions unoccupied by circles) in the regular hexagon 

a = side of regular hexagon 

R = radius of the smallest circle circumscribing all the circles the most densely packed in the regular hexagon 


L = minimum length of one turn of string (i.e. thin, perfectly flexible, and inextensible) wrapped around all the 
packed circles 


Now, consider the following cases of the value of n i.e. number of identical circles tangent to each side of regular 
hexagon as follows 


Case-1: n = 2 (i.e. two identical circles are tangent to each side of regular hexagon) 


Consider total N number of identical circles each of radius r the most densely packed in a regular hexagon 
ABCDEF of side a such that two identical circles are tangent to each side say side ED (as shown in the figure- 
16(a) below). 


Figure-16: (a) Two identical circles labeled with 1, 2 are tangent to (each) side ED of regular hexagon 
ABCDEF (n=2) (b) smallest circle circumscribing all the packed circles (c) minimum length of one turn of 
flexible & inextensible string wrapped around the identical circles most densely packed in a regular hexagon. 


In this case, unknown parameters N, N,,,a,R, and L (as defined above) can be easily determined as given below. 
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In order to find out the number of identical circles and number of voids, the centres of circles lying on the sides 
& vertices of internal regular hexagons are counted in CW circular fashion (as shown in the figure-17) as follows 


N = Total number of identical circles packed in regular hexagon 
=1+6 


It is worth noticing that 6 voids are always present at 6 vertices/corners of the 
outermost regular hexagon ABCDEF (see figure-17). 


N,, = Total number of voids (unoccupied regions) in regular hexagon 
= 6 + 2(6) Figure-17: The centers of 

circles lie on inner hexagon 
F — which can be easily counted in 

Ie right BiGAs We have 2 Abia = CW circular fashion around 


16(a) above) the central circle 


a “ = 30°, C,G =r (see the figure- 


AG r 
tan ZAC,G =—~ > AG=(C,G tanZAC,G =rtan30° = — 


CG V3 
2r 2r 
a = Side of regular hexagon ABCDEF = AB = AG +GH + HB = GH + 2(AG) = 2r + == (2 -1)(2r) + = 
v3 v3 
R = Radius of smallest circumscribing circle = ON = r+2r=r+(2—-1)(2r) (see figure 16(b) above) 
In the figure-16(c) above, We have 2XC,S = WS =° =" Cs=Gx=r 
=> ArcSX = C,X (ZXC,S) =r(=) ae 
a, ee N67 6 
Yr Tr 
XY =arcXS +ST +arcTY =ST + 2(arc XS) = 2r +2 (=) =r +> (see figure 16(c) above) 


mr mr 
L = Minimum length of one turn of string = 6(XY) = 6 (2r + =) =6 (c —1)(2r) + 5) 


Case-2: n = 3 (i.e. three identical circles are tangent to each side of regular hexagon) 


Similarly, consider total N number of identical circles each of radius r the most densely packed in a regular 
hexagon ABCDEF of side a such that three identical circles are tangent to each side say side ED (as shown in the 
figure-18(a) below). 


(b) (c) 


Figure-18: (a) Three identical circles labeled with 1, 2, 3 are tangent to (each) side ED of regular hexagon 
ABCDEF (n=3) (b) smallest circle circumscribing all the packed circles (c) minimum length of one turn of 
flexible & inextensible string wrapped around the identical circles most densely packed in a regular hexagon. 


In this case, the unknown parameters N, N,,a,R, and L can be easily determined as given below. 
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In order to find out the number of identical circles and number of voids, the centres of circles lying on the sides 
& vertices of internal regular hexagons are counted in CW circular fashion (as 
shown in the figure-19) as follows 


N = Total number of identical circles packed in regular hexagon 
=1+6+4+12 


It is worth noticing that 6 voids are always present at 6 vertices/corners of the 
outermost regular hexagon ABCDEF (see figure-19). 


N,, = Total number of voids (unoccupied regions) in regular hexagon Figure-19: The centers of circles 
=~ 6426412 lie on inner hexagons which can 
=6+2(6+12) be easily counted in CW circular 


fashion around the central circle. 


an a = 30°, C,G =r (see the figure- 


In right AC,GA, We have ZAC,G = 
18(a) above) 


AG T 
tan ZAC,G =—~ > AG =C,G tanZAC,G =rtan30° = — 
1 C1G 1 1 af 


2r 2r 
a = Side of regular hexagon ABCDEF = AB = AG + GH + HB = GH + 2(AG) = 2(2r) + a = (3 -1)(2r) + a 


R = Radius of smallest circumscribing circle = ON = r+ 2(2r) =r + (3—-1)(2r) (see figure 18(b) above) 
In the figure-18(c) above, We have 2XC,S =“W2 =° =" Cs=Gxe=r 
=> ArcSX = CX (ZXC,S) =r(Z) ue 
a Pn NGI iG 
Tr mr 
«XY =arcXS +ST +arcTY = ST + 2(arc XS) = 4r +2 (=) = 2(2r) +> (see figure 18(c) above) 


mr mr 
L = Minimum length of one turn of string = 6(XY) = 6 (2(2r) + =) =6 (c —1)(2r) + 5) 


Case-3: n = 4 (i.e. four identical circles are tangent to each side of regular hexagon) 


Similarly, consider total N number of identical circles each of radius r the most densely packed in a regular 
hexagon ABCDEF of side a such that four identical circles are tangent to each side say side ED (as shown in the 
figure-20(a) below). 


Figure-20: (a) Four identical circles labeled with 1, 2, 3, 4 are tangent to (each) side ED of regular hexagon ABCDEF (n=4) 
(b) smallest circle circumscribing all the packed circles (c) minimum length of one turn of flexible & inextensible string 
wrapped around the identical circles most densely packed in a regular hexagon. 
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In this case, the unknown parameters N, N,,, a, R, and L can be easily determined as given below. 


In order to find out the number of identical circles and number of voids, the centres of circles lying on the sides 
& vertices of internal regular hexagons are counted in CW circular fashion 
(as shown in the figure-21) as follows 


N = Total number of identical circles packed in regular hexagon 
=1+6+12+18 


It is worth noticing that 6 voids are always present at 6 vertices/corners of the 
outermost regular hexagon ABCDEF (see the figure-21). 


N,, = Total number of voids (unoccupied regions) in regular hexagon 
=6+2(6+12+4 18) 


Figure-21: The centers of circles lie 
on inner hexagons which can be 


In right AC,GA, We have 4AC,G =—— = —= 30°, C,G =r (see the asily counted in CW circular 


ZAOB _ 6 


fashion around the central circle. 


figure-20(a) above) 


— Ae = = o= 7 
tan ZAC,G = eG => AG=C,G tanZAC,G =rtan30° = a 
2r 2r 
a = Side of regular hexagon ABCDEF = AB = AG + GH + HB = GH + 2(AG) = 3(2r) + a = (4-1)(2r) + a 
R = Radius of smallest circumscribing circle = ON = r+ 3(2r) =r+(4—-1)(2r) (see figure 20(b) above) 
In the figure-20(c) above, We have 2XC,S = “0 = a =2,GS=OX=r 
=> Arc SX = C,X (£XC,S) =r(c) = 
= me N67 6 
cr Tr 
«XY = arc XS + ST +arcTY = ST + 2(arc XS) = 6r +2 (=) = 3(2r) +> (see figure 20(c) above) 


mr mr 
L = Minimum length of one turn of string = 6(XY) = 6 (3(2r) + =) =6 ( —1)(2r) + 5) 


Case-4: n = 5 (i.e. five identical circles are tangent to each side of regular hexagon) 


Similarly, consider total N number of identical circles each of radius r the 
most densely packed in a regular hexagon ABCDEF of side a such that five 
identical circles are tangent to each side say side ED (as shown in the figure- 
23(a) below). 


In this case, the unknown parameters N,N,,a,R,and L can be easily 
determined as given below. 


In order to find out the number of identical circles and number of voids, the 


centres of circles lying on the sides & vertices of internal regular hexagons 


are counted in CW circular fashion (as shown in the figure-22) as follows —_ Figure-22: The centers of circles lie on 
inner hexagons which can be easily 
counted in CW circular fashion around 


N = Total number of identical circles packed in regular hexagon F 
the central circle. 


=1+6+12+18+4 24 
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It is worth noticing that 6 voids are always present at 6 vertices/corners of the outermost regular hexagon ABCDEF 
(see figure-22 above). 


N,, = Total number of voids (unoccupied regions) in regular hexagon 
=6+2(6+12+4+ 18+ 24) 


Figure-23: (a) Five identical circles labeled with 1, 2, 3, 4, 5 are tangent to (each) side ED of regular hexagon ABCDEF (n=5) 
(b) smallest circle circumscribing all the packed circles (c) minimum length of one turn of flexible & inextensible string 
wrapped around the identical circles most densely packed in a regular hexagon. 


In right AC; GA, We have ZAC,G = a = “ = 30°, C,G =r (see the figure-23(a) above) 


ee = = oT 
tan ZAC,G = Ee => AG =(,G tan ZAC,G = rtan30° = a 
2 2 
a = Side of regular hexagon ABCDEF = AB = AG + GH + HB = GH + 2(AG) = 4(2r) + G =(5-1)(2r) + a 
R = Radius of smallest circumscribing circle = ON = r+4(2r) =r+(5—1)(2r) (see figure 23(b) above) 
In the figure-23(c) above, We have 2XC,S = “0 = on =2,GS=OX=r 
=> Arc SX = C,X (ZXC,S) =r(=) a 
= mi" \6/ ~ 6 
Tr Tr 
XY =arcXS +ST +arcTY = ST + 2(arc XS) = 8r +2 (~) = 4(2r) +> (see figure 23(c) above) 


mr mr 
L = Minimum length of one turn of string = 6(XY) = 6 (4(2r) + =) =6 ( —1)(2r) + <] 


Case-5: n = 6 (i.e. six identical circles are tangent to each side of regular hexagon) 


Similarly, consider total N number of identical circles each of radius r the most densely packed in a regular 
hexagon ABCDEF of side a such that six identical circles are tangent to each side say side ED (as shown in the 
figure-25(a) below). 


In this case, the unknown parameters N, N,,a,R, and L can be easily determined as given below. 


In order to find out the number of identical circles and number of voids, the centres of circles lying on the sides 
& vertices of internal regular hexagons are counted in CW circular fashion (as shown in the figure-24 below) as 


follows 
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N = Total number of identical circles packed in regular hexagon 
=1+6+12+18+ 24+30 


It is worth noticing that 6 voids are always present at 6 vertices/corners of 
the outermost regular hexagon ABCDEF (see figure-24). 


N,, = Total number of voids (unoccupied regions) in regular hexagon 
=6+2(6+124+18+ 24 + 30) 


ZAOB _ 6 


In right AC,GA, We have ZAC,G = a" = 30°, C,G =r (see the 
2 


figure-25(a) below) 


AG 
tan 2ACG == = AG = C,G tan ZAC,G = rtan30° = a 
Figure-24: The centers of circles lie on 
inner hexagons which can be easily 
counted in CW circular fashion around 
the central circle. 
en n=6 eels os Re 
COOoge 242846 is ~ 
( Y Y Y Y ce CY YY Y : WAN AX Nf a. A. ah 
Seeese: oe Ory rrr rr ee soe ee 
AYAAAAOO, (| OOOO yy oon 
VOCCCOOCOCOUCOCOOOCOCOIMOCCCC Coen 
AAA RAD: CODE OO EY LYS 
COCCONSSCOO MIO OOO OOO COSI OOOO NOGA 
CLK YY VOSA COO SX 
COP ARAAY  \ CODON XTRA LY 
) Sees eee8 i Y WY Y Wy 
{ YY YY YY y ABMAAAAL) 
AN I A D ) ( Y a Y Y Y\Y “°y 
weheaea COO 4 
—— a ee IAN A. OT OUD wi 


(b) = {c) 

Figure-25; (a) Six identical circles labeled with 1, 2, 3, 4, 5, 6 are tangent to (each) side ED of regular hexagon ABCDEF 
(n=6) (b) smallest circle circumscribing all the packed circles (c) minimum length of one turn of flexible & inextensible 
string wrapped around the identical circles most densely packed in a regular hexagon. 


a = Side of regular hexagon ABCDEF = AB = AG + GH + HB = GH + 2(AG) = 5(2r) + = (6—-1)(2r) + = 
R = Radius of smallest circumscribing circle = ON = r+5(2r) =r + (6—1)(2r) (see figure 25(b) above) 
In the figure-25(c) above, We have 2XC,S =“ = “ =",¢S8=GX=r 
8 mr 
=> Arc SX = C,X (ZXC,S) =r (=) == 


mr mr 
XY = arcXS + ST +arcTY = ST + 2(arc XS) = 10r +2 (=) =5(2r) +5 (see figure 25(c) above) 
m1 Tr 
L = Minimum length of one turn of string = 6(XY) = 6 (s(2r) + =) =6 ( —1)(2r) + 5) 
Similarly, consider n number of identical circles touching each side say side AB of regular hexagon ABCDEF 


then following the symmetry in the above five cases (i.e. n = 2,3,4,5,6,...) (see the figure-26 below), the 
unknown parameters N, N,, a, R, and L can be easily determined as follows 
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Total n number of identical circles are tangent to each side of regular hexagon ABCDEF 


w Naf \ 
Go00 OQQ----- BAA 


B 
t€ a = 


Figure-26: n number of identical circles labeled with 1, 2, 3,...k, k+1, .... n are tangent to (each) side AB of regular hexagon 


ABCDEF (n=n). 


Total number of identical circles (N) packed in regular hexagon ABCDEF 
N=1+64124+184 244 304. .....0...0....2..Upto n terms 


=146(1 4243444 St... cece eee Upto (nm — 1) terms) 


-146(S-20—*8 


) =a anc = 1) = 3n?— 344 


. N=3n?—-3n+1 echoat eects 20) 


Total number of voids (unoccupied regions) (N,,) in regular hexagon 


N, =64+2(64+124+ 184 244 30+.............upto (n — 1) terms) 
=6+12(114+24+34+4454..............upto (n — 1) terms) 
n—-1)(n-1+1 
= 6+12( er *) = 6 + 6n(n = 1) = 6(n? =n +1) 
. N, = 6(n?—n+1) i tiase ate Ret OL) 


The side (a) of the regular hexagon 


=(n-N@n+2=2 ( +-1) 
a=(n Tr 3 rin 3 


. a= 2r(n+—=-1) pastvisivnsenza(2Z2) 


The radius (R) of smallest circle circumscribing all the identical circles packed in the regular hexagon 
R=r+(n-1)(2r) = (2n-1)r 
« R=(2n-1)r Hisietneeen(23) 
The minimum length (L) of the string 


p= 6(cn- 1)(2r) +3] 96n bale 


. L=2(6n+n-6)r sadeuschasrabesenae (24) 
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7.1 Packing fraction/density of identical circles most densely packed in a regular hexagon 


Consider total N number of identical circles each of radius r most densely packed in a regular hexagon ABCDEF 
of each side a such that n number of circles are tangent to each side of the hexagon (refer to the figure-25 above). 


The total number (1) of identical circles most densely packed in regular hexagon is given by the above Eq(20) 
N = 3n? —3n+1 


The side (a) of regular hexagon is given by the above Eq(22) 
2 ( + : 1) 
a=2r(n+—- 
v3 


Now, the packing fraction or packing density (p) of identical circles in regular hexagon is given as follows 


_ Area occupied by N number of identical circles each of radius r 


Area of regular hexagon with each side a 


-NG@r?) Grt—3n4¢1)er? —  Gn?=—3n+1)n _ Gn? —3n+1)nv3 
“BV3 5 . a3 1 2 BGR =D 21) 6Ge@e +4) 
2 393 (ar (n+ 25-1)] 


2_ 
= clini as ee © cf 


.—p —_____—_, 
6 (V3(n — 1) +1) 
It is clear from the above Eq(25) that the packing fraction/density of identical circles of a finite radius most densely 


(hexagonally) packed in a regular hexagon depends on the number n of identical circles tangent to each side but 
independent of the radius of identical circles. 


The variation of packing fraction or density (p) w.r.t. to number of most densely packed identical circles (n) 
tangent to each side of a regular hexagon has been shown in the figure-27. 


MATLAB Code to plot variation of packing density (p) : : . = SeNS'D 
w.r.t. number of circles (n) tangent to each side 


clc; 

clear all; 

close all; 

syms xX y 

% x=[0:0.2:6.5]; 

y=(pi*sqrt(3)/6)*((3*x.%2-3*x+1) /(sqrt(3)*(x-1)+1)%2); 
fplot(x,y,[2, 100]); 

title('\rho vs n'); 

xlabel('Number of circles (n) \rightarrow'); 
ylabel('Packing density (\rho)\rightarrow' ); es 
grid on; 


Packing density (p) + 
° 
8 
7 


° 
s 


10 20 30 40 50 60 70 80 90 100 
Number of circles (n) + 


It is clear from the above graph that the packing density (p) 


of identical circles most densely packed in a regular hexagon Figure-27: Variation of packing density (p) w.r.t. number of 


: . ; . : most densely packed identical circles (n) tangent to each side 
increases with the number of identical circles (n) tangent to of regular hexagon. 


each side and reach to a maximum value when n > ©0. 


The minimum value of packing density Pin of identical circles most densely packed in a regular hexagon can be 
determined by substituting n = 2 in the above Eq(25) as follows 
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= 2B (sar s@ +) mB 7 )-8 7(v3 -1) _ 7n(2v3 - 3) 


m6 \wse-n+1)') © \ws+0)') © \(as+003-0) 12 


71 
Spa = Tz (2¥3 - 3) = 0.8505106310376239 rns i 


VP min = 85.05 


It is clear from the above graph that the packing density (p) has its minimum value ~ 0.85 for n = 2. 


7.2 Maximum packing fraction/density of identical circles packed on an infinite plane 


In order to obtain the maximum packing fraction or density (p), the identical circles have to be packed in a regular 
hexagon. An infinite plane can be assumed to be a regular hexagonal plane with side tending to infinity. In this 
case, the number of identical circles touching each side of infinite regular hexagon will tend to infinity i.e. 2 — 00 


Therefore the maximum packing fraction or density of identical circles (0 ,,) of a finite radius packed over an 
infinite plane can be determined by taking limit of packing density p given from Eq(25) as n — © as follows 


3,1 
ae 3n? — 3n+1 )= mv3/ 3 -ataz =“ ( 3-0+0 )- 
1 


Pmax = lim lim 
n-0o 6 (V3(n-1) +1)" no 6 (v3(1-2)+4) 6 (V3 — 0) +0)” 6 
Ttv3 
 Pmax = “aks = 0.9068996821171088 Semananencl lt) 


6 


%Pmax = 90.69 


It is clear from the graph in above figure-27 that the packing fraction/density (p) of most densely packed identical 
circles in a regular hexagon starts from its minimum value ~ 0.85 for n = 2 number of identical circles (n) 
tangent to each side of regular hexagon and reaches its maximum value pmax = a8 = 0.9069 when n > i.e. an 


infinite plane. 


It is worth noticing that pmax = ma = 0.9069 is maximum possible packing fraction/density of identical circles of a 


finite radius packed over an infinite plane. 


Alternatively, the maximum packing fraction/density pmax of identical circles packed 
over an infinite plane can be obtained by considering a circle of finite radius 7 and 
centre O surrounded by 6 external tangent identical circles with centres A, B, C, D, E 
and F packed in the hexagonal pattern (as shown in the figure-28). Consider a regular 
hexagon ABCDEF by joining the centres of 6 identical circles surrounding the 
identical circle with centre O. 


Therefore the maximum packing fraction or density of identical circles (Omax) is 


determined as follows 


Total area occupied by circles in regular hexagon ABCDEF Figure-28: The identical 
Pmax = circles are most densely 
Area of regular hexagon ABCDEF packedlover an iulinite plane, 
_ 3@rr?) oa _ av3 
Pmax 33 a3 6 


(er)? 
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8. 2D packing of circles in the regions bounded by two external tangent circles and the 
circumscribed circle 


When two external tangent circles of radii say a, & b are inscribed in a third circumscribed circle of radius R then 
the unoccupied region is divided into two regions which are always unequal (i.e. minor and major) except when 
R=a+b (See the above figures 2 & 3 or the figure-29 below). 


The circle packing in a plane region, bounded by three circles out of which two are external tangent circles 
touching the third circle internally (as shown in the figure-29 below), can be easily analysed and formulated by 
substituting a = 7; (i.e. variable/input radius) and Cin = 741 (i.e. radius of next packed circle) in the Eq(1) (as 
derived in the above section 1.1) to analytically compute the radii of packed circles using the generalized 
recurrence formula given by as follows 


r,bR ((R — b)r, + bR — 2,/bRr,(R — b — r)) 


Ti41 = V Ti41 < Tj see eee eee eee (28) 


((R + b)r, + BR) — 4r;bR? 


r,bR ((R — b)r; + BR + 2,/bRr{(R — b — r)) 


Tia. = V T; < Ti44 <R-b aac Woarenaresaya (29) 


((R + b)r; + BR) — 47,bR? 


Where, 7, = a is the radius of one of two external tangent circles i.e. initial circle and b is the radius of second 
external tangent circle i.e. reference circle. Both the external tangent circles internally touch the third (i.e. largest) 
circumscribing circle of radius R and divide unoccupied region into two parts i.e. region-1 & region-2 (see figure- 
29 below). 


1) The first recurrence relation i.e. Eq(28) always gives the diminishing or decreasing radii of packed circles in 
minor region-1 (i.e. labelled with grey colour) and therefore Eq(28) must be used for computing radii of circles 
packed in minor region-1 for given values of a, b, and R (as shown in the figure-29). 


Reference circle 


Figure-29: Circle packing in the plane regions 1 and 2 (i.e. 
labeled with grey and yellow colors) bounded by three 
mutually tangent circles with radii a, b and R 
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2) The second recurrence relation i.e. Eq(29) gives the increasing radii up to a certain maximum value such that 
Yi44 <R — b of packed circles in major region-2 (i.e. labelled with yellow colour in the figure-29 above) and after 
attaining maximum value, the radius 7;,, starts oscillating from 7; to its maximum value recursively. Therefore 
the Eq(29) must be used for packing circles of only increasing radii in major region-2 such that 7, < 74, < R—b 
and after 7;4, reaches its maximum value, Eq(28) must be used to further pack required number of circles (of 
decreasing radii) in the major region-2. 


It is worth noticing that in circle packing into major region-2 using the Eq(29), 


AbDR(R — b) 


Raby 


> N41 = 

The above case shows that the radii 7, and 7;,, of two circles consecutively packed into major region-2 (i.e. 
labelled with yellow colour) will be equal and both will be tangent to the diameter line of circumscribing circle 
of radius R (as shown in the figure-30 below). In this case, the maximum possible value of 7,1 is given as 


_ ABR(R — b) _ 4bR(R — b) 
(Ti+1) max = ~(R+b)2_ i= “(R+b)2_ 


In this case of circles packing into major region-2, after attaining 


; ; 4bR(R-b) 
the maximum value i.e. (1341) max = (R+b)? 


will not oscillate because 7;4, = 7; i.e. radii of two consecutively 
packed circles become equal. This is Non Resonance (NR) 


, the radius 744 


4bR(R — b) } 


condition of the above recursive Eq(29). t= Rape / 


While in all other cases of circle packing in major region-2 using 
recursive Eq(29), the radius 7:4, of packed circle will always 
4bR(R-b) 
(R+b)2 
consecutively packed circles are distinct i.e. 714, #7; . This is 


oscillate from 7; to (741) max # because the radii of two 


Resonance (R) condition of the above recursive Eq(29). 


i) The value of 7:4; obtained from the above recursive Eq(29) 


: : : : 4bR(R-b : ; : : 
will not oscillate if and only if 7; = sca Figure-30: The radius 1;,, of circle packed into 
(R+b) major region-2 (i.e. labelled with yellow color) will 
= . : not oscillate if and only ifr; = ue : 
ii) The value of 7;4, obtained from the above recursive Eq(29) (Rea) 
4DR(R-b) 


will always oscillate if 7; # ipebe 
3) In the above recurrence formula i.e. Eq(28) and Eq(29), one of two radii a and b is changed to r, (while other 
is kept constant/fixed) and 7;,, takes the value of radius of a newly packed circle 7; recursively until the desired 
number of circles are packed in the bounded plane region-1 or region-2 or both (see figure-29 above). 


4) In order to pack n, number of circles into minor region-1 and n, number of circles into major region-2 of area 
unoccupied by two external tangent circles of given radii a, and b which internally touch the third circle of radius 
R, the following MATLAB code can be used to determine the radii of all the desired number of circles packed 
into two unoccupied regions 


cle; 

clear all; 

close all; 

R=input('Enter the radius of circumscribing circle (in required unit) = '); 
a=input('Enter the radius of variant external circle. a = '); 
b=input('Enter the radius of reference external circle, b = '); 
n1=input('Enter the number of circles to be packed in minor region= '); 
n2=input('Enter the number of circles to be packed in major region= '); 
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disp('(1) Radii of circles packed into minor region'); 
x=[]3 
for i=1:n1 
r(1)=a; 
r(i+1)=double((r(i)*b*R*(r(i)*(R-b)+b*R-2*sqrt(r(i) *b*R*(R-b-r(i)))))/((r(i)*(R+b) +b*R)%2-4*r (i) *b*R2)) 5 
A=[x, real(r(i+1))]; 
disp(A); 
end 
disp('(2) Radii of circles packed into major region'); 
for i=1:n2+10 
s(1)=a; 
s(i+1)=double((s(i)*b*R*(s(i)*(R-b)+b*R+2*sqrt(s(i)*b*R*(R-b-s(i)))))/((s(i)*(R+b) +b*R)*%2-4*s (i) *b*R%2)) 5 
M=[x,s(i)]; 
mi=max(M) ; 


end 

for i=1:n2+11 
s(1)=a; 
s(i+1)=double((s(i)*b*R*(s(i)*(R-b)+b*R+2*sqrt(s(i)*b*R*(R-b-s(i)))))/((s(i)*(R+b)+b*R)*2-4*s5 (i) *b*R%2)) 5; 
M=[x,s(i)]; 
m2=max(M) ; 

end 

m=max(m1,m2) ; 


for i=1:n2 
Try ipalalet essen oan eee Ze an iy ea 
if t(i+1)>t(i)+10*(-17) && t(i+1)<m+10(-17) 
Be [x, real(t(i#1))]5 
disp(B); 


end 

end 
c_max=doub1le((m*b*R*(m*(R-b)+b*R+2* sqrt (m*b*R* (R-b-m))))/((m*(R+b)+b*R)*2-4*m*b*R42) ) 5 
c_min=doub1le((m*b*R*(m*(R-b)+b*R-2*sqrt(m*b*R*(R-b-m))))/((m*(R+b)+b*R)*2-4*m*b*R‘2) ) ; 


disp(m) ; 
disp('Set-1'); 


for i=1:n2-n 
u(1)=c_max; 
u(i+1)=double((u(i)*b*R*(u(i)*(R-b)+b*R-2*sqrt(u(i)*b*R*(R-b-u(i)))))/((u(i)*(R+b)+b*R)*2-4*u(i) *b*R4%2) ) 5 
C=[x,real(u(i))]; 
disp(C); 
end 
disp('Set-2'); 


for i=1:n2-n 
v(1)=c_min; 
v(i+1)=double((v(i)*b*R*(v(i) *(R-b)+b*R-2*sqrt(v(i)*b*R*(R-b-v(i)))))/((v(i)*(R+b)+b*R)*2-4*v(i) *b*R4%2) ); 
D=[x,real(v(i))]; 


disp(D); 
end 


Note: In the above MATLAB code, two sets i.e. Set-1 and Set-2 of radii of circles packed into the major region 
are obtained out of which only one set of non-repetitive radii is taken into consideration while the other set gives 
the same values of radii of circles as packed into minor region. Hence care is taken to select correct set of values. 


Conclusions 


In this paper, the problems of circles packed into the region bounded by the circular arcs, and the straight lines 
including square, semi and quarter circles as well as packing of identical circles on infinite plane have been 
explained and solved in details. The generalized formula have been derived and established for analysing and 
formulating the circle packing in 3D space as well as on the spherical surface. The simulated values have been 
compared and verified by analytic formula. The mathematical relations, derived here, can be applied for solving 
the problems of dense packing of circles and spheres in 2D & 3D containers in industrial applications, and finding 
the planar density on crystallographic plane. 
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